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1. Introduction 

The TTTT scattering amplitude is one of the fundamental observables in low energy 
particle physics and has been since long the subject of many studies. The pions 
are the lightest strongly interacting particles and have thus a special status. Their 
properties are also strongly influenced by the chiral symmetry present in the limit of 
massless quarks in Quantum Chromo Dynamics (QCD). 

In this paper we add one more step in the discussion of pion properties. We 
calculate tttt scattering in three flavour Chiral Perturbation Theory (ChPT) in order 
to complete the connection between the Kj^^ form factors and tttt scattering following 
from chiral symmetry. We stay in the isospin limit here and neglect electromagnetic 
effects. But first we give a short historical overview of chiral symmetry relevant to 
TTTT scattering. 
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Chiral symmetry was introduced a long time ago and used in the form of current 
algebra and the PCAC assumptions (partial conservation of axial currents.) Wein- 
berg [1] used these methods to derive a result for tttt scattering valid to lowest order 
in meson masses and momenta. It was later realized that the assumptions of ana- 
lyticity used in many PCAC type of analyses were not always true. This allowed to 
calculate often the leading nonanalytic corrections to the lowest order PCAC results. 
This line of work within the PCAC methods has been reviewed in Ref. [2] where 
also references to earlier work can be found. The more modern method of using 
chiral symmetry in the form of an effective field theory was introduced by Weinberg 
[3] and systematized by Gasscr and Leutwyler [4]. This is now known as ChPT. 
They performed the full next-to-leading order (NLO) calculation of tttt scattering 
[5] as the first major application. The parameters necessary for determining the 
TTTT scattering lengths at threshold had to be taken from the D wave experimental 
results. A reasonable agreement with the experimental result was obtained. Gasser 
and Leutwyler also extended ChPT to the three flavour sector including the strange 
quark in addition to the up and down quark in the ChPT formalism [6]. This allowed 
to determine the parameters necessary for prediction tttt scattering to be determined 
from the absolute values of the form factors in Kf^ decays. Note that the phase and 
the absolute values are in principle separately measurable quantities there such that 
the relation between the phase as determined by the relevant tttt scattering phase 
and the absolute value is not a trivial prediction. These NLO calculations of K^^ de- 
cays were performed and again led to a reasonable agreement with the known values 
[7, 8] . The three flavour expression for tttt scattering was first calculated in [9] and 
later independently in [10]. 

The first step at next-to-ncxt-to-leading order (NNLO) was done in Ref. [10] 
where the dispersive part of the amplitude was determined. The full two loop calcu- 
lation in two fiavour ChPT was performed in Ref. [11, 12]. The NNLO calculation of 

has also been performed [13, 14] and was used to study tttt scattering via the two 
fiavour NNLO calculation using the relation between the NLO order parameters in 
two and three fiavour ChPT derived in Ref. [4] . This still leaves an uncertainty since 
in order to have full control at NNLO the corrections to those relations need also to 
be determined. This can be done in principle by integrating out the kaons and eta 
degrees of freedom out of the three fiavour ChPT NNLO generating functional but 
this has not been done so far. Alternatively one can directly calculate the NNLO 
amplitude for tttt scattering in three fiavour ChPT. This is what has been done in 
this paper. 

A different way to describe theoretically tttt scattering is to use the constraints 
from analyticity and unitarity. These lead to many different sum rules but especially 
after crossing properties have been included the resulting set of equations becomes 
very constraining. These are known as the Roy equations [15]. They were analyzed 
extensively in the seventies, see e.g. [16]. The available high energy data together 
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with the Roy equations allowed to describe tttt scattering in terms of two parameters 
usually chosen to be the scattering lengths in the S wave channel for isospin and 
2, Oq and Oq respectively. Gq was then essentially fixed by measuring the difference in 
phase of the K^.^ form factors. This measurement has been dominated for a long time 
by the experiment of Ref. [17]. The result was in disagreement with the Weinberg 
prediction [1] but in borderline agreement with the NLO order one [5] . However, the 
central values was rather different from the ChPT prediction. If this central value 
turned out to be correct, the consequences for the low energy structure of QCD would 
be quite strong. In particular the chiral symmetry breaking might not be driven by 
the simple quark anti-quark condensate, see [18] and references therein. This was the 
origin for the renewed interest in tttt scattering. The analysis of the Roy equations 
has been updated in Ref. [19]. This was then combined with the constraints from 
chiral symmetry in Ref. [20, 21]. The constraints used were the solutions of the Roy 
equations of [19], the ChPT two flavour calculations of tttt scattering [11, 12] and 
the pion scalar form factor [22]. It also used the method of determination of low 
energy constants from sum rules over tttt phases from [23] . This analysis led to very 
constrained predictions for the scattering lengths^. These were confirmed nicely by 
the E865 experiment at BNL [26, 27]. A similar analysis but without the constraint 
from the pion scalar radius can be found in [28]. 

In this paper we will compare the predictions for the tttt scattering lengths in 
three flavour ChPT with the experimental input from decays to NNLO. We flnd 
an acceptable agreement as discussed in Sect. 6. 

In two flavour ChPT it is clear now that spontaneous chiral symmetry breaking 
goes through the quark anti-quark condensate [20, 21]. There is still a possibility that 
the behaviour when also the strange quark mass goes through zero, is qualitatively 
different. This is discussed in the recent work by Stern and collaborators [29, 30]. 
The argument is that large disconnected loop contributions from strange quarks, via 
kaons and etas, can be large, making a convergent three flavour ChPT difficult to 
achieve in the usual sense [30, 31]. This has been studied in some detail in [14] and 
also in the context of the three flavour ChPT calculations of the various scalar form 
factors [32]. The masses and decay constants, see Ref. [33] and references therein, 
showed the possibihty of this behaviour. The various vector form factors calculated 
did not seem to have problems with convergence [34, 35, 36, 37]. It turns out that 
this calculation does not provide much more information than the scalar form factors 
[32] did. Work is in progress to extend the nK scattering also to NNLO. This might 
allow us to shed more light on this issue. 

This paper is organized as follows. In Sect. 2 we give a short overview of ChPT 
and the references for the methods of NNLO calculations. In Sect. 3 the general 

^The conclusions of that analysis have been challenged in [24], the reply of the authors of [20, 21] 
can be found in [25]. 
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properties of the tttt scattering amplitude are described and the quantities used later 
defined. Sect. 4 gives an overview of our main result, the calculation of the tttt 
scattering amplitude to NNLO in three flavour ChPT. We also present here some 
plots showing the importance of the various contributions. The inputs we use to 
do the numerical analysis are described in Sect. 5. The main numerical analysis is 
presented in Sect. 6 and we give our main conclusions in Sect. 7. The appendices are 
devoted to giving references about the various loop integrals used in this work and 
the exphcit expressions at NNLO for tttt scattering. 



2. Chiral Perturbation Theory 



ChPT is the effective field theory for QCD at low energies introduced by Weinberg, 
Gasser and Leutwyler [3, 4, 6]. Introductory lectures can be found in Ref. [38]. This 
leads to an expansion in quark masses and meson momenta generically labeled p and 
assumes ~ since for an on-shell meson = m^. The ChPT formahsm exists 
both for two light flavours, up and down, referred to as SU (2) ChPT, and for three 
light flavours, up, down and strange, referred to as SU{3) ChPT. The Lagrangian 
for the strong and semi leptonic mesonic sector to NNLO can be written as 



jC — jC2-\- jC,4 + jCq , 

where the subscript refers to the chiral order. The lowest order Lagrangian is 



>C2 = ^(«X + x+)- 



[2.1) 



(2.2) 



The mesonic flelds enter via 
f iM \ 



u — exp 



F0V2, 



M 



(2.3) 



and the quantity also contains the external vector (i>^) and axial-vector (a^) 
currents 

u^j, = i{u^d^j,u- d^j,uu^ -iu^r^u + iul^u^) , lfx.{rij) = v^-{+)a^j,. (2.4) 
The scalar (s) and pseudo scalar (p) currents are contained in 

X± = u^xu^ ± ux^u , X = (s + ip) . (2.5) 
The or NLO Lagrangian, £4, was introduced in Ref. [6] and reads 
£4 = Li{uXf + L2{u^u,){uV) 
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+L-j{x-? + U{X+X-) 

+H,{F^,F^^^ + F/^^Ff^-^^) + H,{xl - xi)/4 ■ (2.6) 

The Lg and Lio terms introduce also the field strength tensor 

Fi.^""^ = dAr). - dJir), - t [/(r)„ /(r).] . (2.7) 

The two terms proportional to Hi and H2 are high energy contact terms and are not 
involved in physical amplitudes, Lg and Liq play only a minor role for the quantities 
discussed in this paper. 

We quote the schematic form of the NNLO Lagrangian in the three fiavour case 

£6= E C^0^ (2.8) 
1=1,94 

and refer to [39] for their explicit expressions. The last four terms are contact 
terms [39]. 

The ultra-violet divergences produced by loop diagrams of order and cancel 
in the process of renormalization with the divergences extracted from the low energy 
constants Lj's and Ci's. We use dimensional rcgularization and the modified minimal 
subtraction (MS) version usually used in ChPT. An extensive description of the 
regularization and renormalization procedure including the freedom involved can be 
found in Refs. [12] and [40]. 

The subtraction of divergences is done explicitly by 

= M'-'itA + LKi,)] (2.9) 

and 



'"''^ c[(/.) - (rS^) + t\'\i,)) a - rf A^] (2.10) 



F2 

where c and A are defined by 

lnc = -i[ln47r + r'(l) + l] , (2.11) 

A= -. (2.12) 

The coefficients f j, ^ and ^ are constants while the Pj^^ 's are hnear combinations 
of the LI{ijl) 's. Their explicit expressions can all be found in [40] where they have 
been calculated in general. The NLO divergences were first calculated in Ref. [4, 6] 
and the doubles poles at NNLO first in Ref. [41]. 
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3. The TTTT amplitude: general properties 

The TTTT scattering amplitude in all the relevant channels can be written as a function 
A{s,t,u) which is symmetric in the last two arguments: 

A {7r%p,)AP2) ^ 7r^(P3)7r''(p4)) = 5'''5'^A{s, t, u)+5'''5'''A{t, u, s)+rt''M(ii, t, s) . 

(3.1) 

s,t,u are the usual Mandelstam variables 

s = [pi + t = {pi- psf and u = {pi- p^f . (3.2) 

The various isospin amplitudes can be written in terms of this function as 

T°(s, t) = 3A(s, t, u) + A{t, u, s) + A{u, s, t) , 
r^(s, t) = A{t, u, s) - A{u, s, t) , 

T'^{s,t) ^ A{t,u,s)+A{u,s,t), (3.3) 
where the kinematical variables t, u can be expressed in terms of s and cos 9 as 

t^-]-{s-Aml){l-cose), u^-hs-Aml){l + cose). (3.4) 
The various amplitudes can be expanded in partial waves via 

oo 

T\s, t) = 327r Y,{21 + l)P^(cos ^)4(s) . (3.5) 

1=0 

Near threshold these are expanded in terms of the threshold parameters 

4 = q'' [al + hW + 0{q')) , q' = \{s- Ami) ■ (3-6) 
Below the inelastic threshold the partial waves satisfy 



Imti(s) = a(s) 



\ a{s) = Jl-^. (3.7) 



In this regime the partial waves can be written in terms of the phase-shifts as 

tlis) = I ^ {e'-''^^'^ - 1} . (3.8) 

In ChPT the inelasticity only starts at order p^. The p^ result has only nonzero 
items for ^q, t\ and ^q. As a consequence the imaginary parts for all other partial 
waves starts only at order p^. In [18] it has been shown that thus up to order p^ the 
amplitude can be written as 

A{s,t,u) = C{s,t,u) + 327r(^^V'{s) + ^{s- u)V\t) + ^{s - t)V\u) 

+ lv\t) + lv\u) - W\s)\ . (3.9) 
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The function V\s) have a polynomial ambiguity but we have chosen to keep it in this 
more general form. The V^{s) contain the singularities in the nn amplitudes from 
intermediate states with isospin I in the various channels. They obey the relations 

Imy°'2(s) = Imf/{s) , lmV\s) = _\ ^ lmtl{s) . (3.10) 

The polynomial C{s, t, u) can be written using s + t + u — 4m^ in the form 

C{s, t, u) = ci + C2S + c^s^ + C4(t — m)^ + 056'^ + CQs{t — uY . (3.11) 



4. ChPT results 

4.1 Two Flavour ChPT 

The lowest order was derived by Weinberg [1] a long time ago and corresponds to 

= (4-1) 

in (3.11). The next order was derived by Gasser and Leutwyler [5] and the full 
calculation of order was performed in [11, 12]. 

4.2 Three Flavour ChPT 

The lowest order is identical to the two-flavour case of Eq. (4.1). The order nn 
scattering amplitude in three flavour ChPT was first published in [9] App. A. It can 
also be found in [10]. Notice that the contribution from and L3 is missing in [9]. 
Our result at this order is in full agreement with the corrected version. We have 
expressed the result in terms of the functions defined in Eq. (3.9) in App. B. 

The p^ expressions we present arc those corresponding to the p"^ results expressed 
in the physical masses and decay constants. The order p^ expression is our main 
result. Expressed in the polynomial C{s,t,u) and the functions y°'^'^(s) it is shown 
in App. B. 

4.3 A first numerical look 

In this subsection we present a first look at the numerical results for the two loop 
amplitudes. We choose as input the pion decay constant, the charged pion mass, an 
averaged kaon mass with electromagnetic effects removed and the physical eta mass. 

= 92 A MeV , = = 139.56995 MeV , 
rriK = 494.53 MeV , = 547.3 MeV . (4.2) 

The subtraction scale /i — 770 MeV is used throughout the paper unless otherwise 
mentioned explicitly. 
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We first compare the SU (2) and SU (3) ChPT results for tlie tliree main partial 
waves for all low energy constants set to zero at the scale 770 MeV. So we set 
L\ = CI = II = = Q and compare the pure loop results for the SU (2) and SU (3) 
ChPT. 

loops only Lj' = Cj' = 




-0.05 I ' ' ' ' ' ' ^ 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 

s [GeV^] 



Figure 1: Results for the t[] partial wave: The SU{2) ChPT results of [12] for all LECs 
equal to zero, labeled 2F, compared to the results for the SU{'i) ChPT results also for all 
LECs equal to zero at ;U = 770 MeV, labeled 3F. 

In Fig. 1 we have plotted the partial wave amplitude t^{s) Notice the extremely 
small difference between the two results, showing that this channel is obviously dom- 
inated by the pion loops and the kaon and eta have only a fairly small effect. 

For the t\ partial wave, shown in Fig. 2 the order results are very similar in 
both cases again but there is somewhat more difference in the order contributions, 
this channel has a much stronger effect from the LECs, see below, so this difference 
does not play much of a role. 

The amplitude, shown in Fig. 3 has a somewhat more surprising difference. 
Again the order results for the SU (2) and S'f/(3) cases are very similar but the 
pure loop correction is rather different. We will come back to this later. 

That the difference at order p^ would be small was of course expected. In [12] 
the contribution from the p^ loops from kaons and etas was estimated and found to 
be very small. 

5. Solution of the Roy Equations and Other Inputs 

For the experimental values of the tttt scattering phase-shifts we use the results from 
the extensive analysis of the Roy equations done in [19]. This analysis has been 
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loops only ^ = C\ = 

0.04 

0.03 

0.02 

0.01 
rj- 
-0.01 
-0.02 
-0.03 
-0.04 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 

s [GeV^] 

Figure 2: Results for the t\ partial wave: The SU{2) ChPT results of [12] for ah LECs 
equal to zero, labeled 2F, compared to the results for the SU (3) ChPT results also for all 
LECs equal to zero at ;U = 770 MeV, labeled 3F. 

loops only ^ = C\ = 

0.06 
0.04 
0.02 


-0.02 
-0.04 
-0.06 
-0.08 
-0.1 
-0.12 
-0.14 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 

s [GeV^] 

Figure 3: Results for the 4 partial wave: The SU{2) ChPT results of [12] for all LECs 
equal to zero, labeled 2F, compared to the results for the SU{3) ChPT results also for all 
LECs equal to zero at ^ = 770 MeV, labeled 3F. 

checked with somewhat relaxed input assumptions in [28]. These equations and the 
input are constrained in such a way that the tttt phase-shifts are determined as a 
function of two input parameters chosen to be the scattering lengths and Og. 
The solution of the Roy equations has been used in [20, 21] together with the 
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constraints from the SU (2) ChPT expression to obtain a precise prediction of Oq 
and Cq. The input used there was the pion scalar radius evaluated using dispersive 
methods and the estimates of the low energy constants of [12] , together with the 
SU{2) order p^ calculations of tttt scattering [11, 12] and the pion scalar radius [22]. 

The conclusions of that analysis have been challenged in [24], the reply of the 
authors of [20, 21] can be found in [25]. 

For our SU{3) ChPT results we use as inputs the masses and decay constants 
given in Sect. 4.3 and a subtraction constant /j, — 770 MeV. We work in the limit of 
exact isospin. 

It should be mentioned that these numbers were obtained using a general refit 
of all the other LI with a fixed and Lg as input. The inputs were the absolute 
values of the Ki4 form factors as measured by the E865 experiment [26, 27], decay 
constants and meson masses. The fits correspond precisely to those of fit 10 in [42] 
but with different input values for and Lg. 

For the constants at order p^ we need to use input values using various estimates. 
The method used was proposed at NLO in Ref. [43, 44] and references therein. 
We do not estimate here the NLO constants this way but only NNLO. The places 
where comparisons with experiment are available are in general reasonable agreement 
with these estimates. The estimates are obtained by including the main resonance 
exchange contributions and putting the p^ part of these amplitudes equal to the 
contribution from the C[. This procedure is obviously subtraction point dependent 
and is normally only performed to leading order in the expansion in l/Nc, with Nc the 
number of colours. Alternative approaches exist but we will not discuss them here. 
Recent papers addressing this type of issues are [45. 46, 47] and references therein. 
A systematic study of this issue is clearly important, for our present purpose the 
estimates seem sufficient. 

The estimates from resonance exchange for the masses and decay constants are 
the most uncertain. These are discussed in [33] and [45]. For the numerical results 
used here they have been put to zero, the naive size estimate of [33] led to extremely 
large NNLO corrections. The estimates of the Ke4 amphtudes can be found in [14] 
after the work of [48]. The effect of varying these was studied in [14] and found to 
be reasonable. 

The estimates of the p^ contributions to tttt scattering we use are those of [12]. 
These lead to the contributions to the various threshold parameters given in Table 1. 
The uncertainty on these is quite considerable but probably within a factor of two, 
this is also discussed in Ref. [21]. Similar resonance estimates of nn scattering can 
be found in [49]. 

In order to be able to perform a study of the dependence on and Lg we 
have redone the fits to the masses, decay constants and the Ki^ absolute values of 
the form factors with a range of values for and Lg. This is similar to the part 
discussed in [14] but now with the newer experimental input [26, 27] included. The 
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p2 


p\ LI = 


p6, L\ = CI = 


CI only 


Ref. [21] 


Uq 


0.159 


0.041 


0.011 


0.001 


0.220 ±0.005 


W 


0.182 


0.075 


0.016 


0.004 


0.276 ± 0.006 


10 al 


-0.454 


0.037 


0.015 


-0.004 


-0.444 ±0.010 


10 6^ 


-0.908 


0.111 


0.029 


-0.011 


-0.803 ± 0.012 


10 a} 


0.303 


0.022 


0.025 


0.000 


0.379 ± 0.005 


10 6^ 


— 


0.005 


0.033 


0.000 


0.057 ±0.001 


102 aO 


— 


0.121 


0.050 


0.001 


0.175 ±0.003 


102 feo 


— 


-0.040 


-0.005 


0.005 


-0.036 ±0.002 


lO^ai 




0.492 


-0.187 


-0.003 


0.170 ±0.013 


10^62 




-0.234 


-0.136 


-0.045 


-0.326 ±0.012 


10^4 




0.20 


0.27 


0.15 


0.560 ±0.019 


lOHl 




-0.15 


-0.22 




-0.402 ±0.018 



Table 1: The various contributions to the threshold parameters defined in (3.6) with the 
L[ = at /U = 770 MeV. The CJ" contribution is shown separately using the estimates of 
[12]. The threshold parameters are given in the corresponding power of m^.,.. Note that 
we have given a\ and b\ always with the same power of ten. 



fits correspond exactly to the fit labeled fit 10 in [42] but with various values of L\ 
and Lg as input. These were already used in Ref. [32] to compare with the scalar 
form factors. There a general preference was found for the region L'^ — 0.0003. 

In that region the corrections to the pion scalar form factor at zero were fairly small 
as well as a good agreement with the pion scalar radius was obtained. It should be 
kept in mind that all the other are varied together with L\ and Lg in order to fit 
the mentioned quantities. Reasonable fits were obtained for most values of L\ and 
Lg. Varying L4 and Lg without the correlated changes in the other L\ would lead to 
much larger variations than the ones shown below. 

6. Numerical Analysis 

In order to check convergence let us first look at the various contributions with all 
the low energy constants set to zero at a scale = 770 McV. These are shown in 
Table 1. The angular integrals have been performed using both a 5 point and a 8 
point Gaussian integration over cos 9. In addition the fits were performed numerically 
over a range of above threshold. The numerical errors on the slopes h\ for £ = 2,3 
are of the order of the last digit shown. For all others this error is below the accuracy 
given. For comparison we have also given the results of Ref. [21] in the last column. 

We will now compare with the full analysis of tttt scattering performed with the 
use of the Roy equations and the SU{2) ChPT results of [20, 21]. In principle we 
could redo this work with the SU (3) ChPT results as constraints instead. We have 
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(a) 



(b) 



Figure 4: The scattering lengths Oq and Oq as a function of the input values L4 and 
Lg with the other L[ simultaneously refitted to the Kf>4 form factors, (a) Oq calculated to 
order and p^. The central value of [21] of 0.220 is also shown, (b) Oq at order and p^. 
The two horizontal planes indicate the allowed region obtained in [21]. 

chosen not to do so, postponing a possible more detailed comparison till after the 
inclusion of ttK scattering results. Since the SU{3) ChPT results contain the pion 
loops which are the main effects of the SU{2) ChPT calculations and the SU{3) 
expressions must reduce to the SU{2) expressions in the limit of a large kaon and eta 
mass and thus satisfy the SU{2) chiral constraints we do not expect such an analysis 
to differ substantially from the one performed in Ref. [20, 21]. 

It can already be seen from Table 1 that the lowest order result together with 
the pure loop contributions only already give quite a good description of the various 
scattering lengths. The effects of including nonzero values for the low energy con- 
stants should explain the difference. Notice that for almost all cases the estimated 
contributions from the constants is rather small. The main effect is thus from 
the p'^ constants L[. We will now study the effects of these when they were fitted to 
other data as described above with fixed values of and Lg as input. 

In Fig. 4(a) we have plotted the result for function of the input values of 

L4 and Lg. The lowest order value. 



is not shown on the plot. The convergence of the series is very good and of similar 
quality as the two flavour result. Taking the result of [21], 



a, 







0.159, 



(6.1) 




a, 







0.220 ±0.005, 



(6.2) 



[21] 
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we see that the agreement is excellent and no new information on and Lg is 
available from this source, this also confirms the prediction for Og. 
The result for is plotted in Fig. 4(b). The lowest order value 



-0.0454 , (6.3) 



is not plotted. The series converges well over most of the L4 and Lg range considered. 
The result should be compared with the result [21], 



= -0.0444 ± 0.0010 , (6.4) 

[zi] 



The two planes in Fig. 4(b) are the error boundaries of Eq. (6.4). We see that in 
order to get agreement we need to go to the front right part of the graph. The point 
with L4 Lg = is outside the error band and has the value 

= -0.0410. (6.5) 

If we take a closer look at the scalar radius results of [32] and especially at Fig. 11(a) 
there, we see that the dispersive and the SU{3) ChPT result for the scalar radius 
are in good agreement at^ 

Ll^Ll- 0.0004 . (6.6) 

The value of the scalar form factor there is in good agreement with the result used 
in [20, 21] of < >f = 0.61 ± 0.04 fm^. Following the line (6.6) in our results leads 
to a virtually constant prediction of 

a^ = -0.0433. (6.7) 

This is in reasonable agreement with the result of [21]. The SU{3) ChPT result thus 
confirms the result of Ref. [21] when the constraint for the scalar radius is taken into 
account. We have not performed a full error analysis for the result (6.7) similar to 
the one performed in Ref. [14], but we expect the errors coming from the various 
uncertainties on the L[ to be similar to the ones quoted there. 

At this point wc have checked the agreement for the two main input parameters 
for the dispersive calculations. How well do the other threshold parameters compare? 
We can show plots similar to the ones for Oq and Oq shown in Fig. 4, but they do not 
provide any essential new information. We have first given in Table 2 the results 
for the various threshold parameters for the input values from fit 10 and fits A,B,C 
as defined in Ref. [32]. This table can be seen as the extension of the one given in 
that reference for the masses, decay constants and scalar radius to the tttt scattering 
threshold parameters. 



^In [32] in addition a small correction to the form factor at zero was required. 
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fit 10 


fit A 


fit B 


fit C 






4 

P 




total 


total 


total 


total 


"o 


0.159 


0.044 


0.016 


0.219 


0.220 


0.220 


0.221 


1.0 

% 


0.182 


0.073 


0.025 


0.279 


0.282 


0.282 


0.282 


10 


-0.454 


0.030 


0.013 


-0.410 


-0.427 


-0.433 


-0.428 


10 62 


-0.908 


0.151 


0.025 


-0.731 


-0.755 


-0.761 


-0.760 


10 a} 


0.303 


0.052 


0.031 


0.385 


0.388 


0.389 


0.389 


10 6} 


— 


0.029 


0.038 


0.067 


0.064 


0.063 


0.063 


102 a" 


— 


0.153 


0.080 


0.233 


0.223 


0.220 


0.221 


102 6° 


— 


-0.040 


0.007 


-0.033 


-0.035 


-0.036 


-0.036 


103 al 




0.327 


-0.106 


0.221 


0.219 


0.218 


0.221 


WHI 




-0.234 


-0.151 


-0.385 


-0.386 


-0.385 


-0.387 


10^ al 




0.20 


0.44 


0.64 


0.62 


0.62 


0.62 






-0.15 


-0.20 


-0.35 


-0.34 


-0.34 


-0.34 



Table 2: The values of the threshold parameters defined in (3.6) for the values of the 
input parameters for the fits 10 [42] and A,B,C [32]. The lowest order values and the 
contributions from the Cf as given in Table 1 are included. The threshold parameters are 
given in the corresponding power of nT?+. Note that we have given a| and 6| always with 
the same power of ten. For fit 10 we also give the three orders separately. 



As can be seen a reasonable agreement is obtained for most threshold parameters 
studied. We have not performed a full error analysis but a first estimate is about half 
the p^ contribution plus the size of the estimated contribution from the C[. Only for 
6o there is a mild discrepancy with these criteria. 

A possible further test can be done by comparing parameters relevant in the sub- 
threshold expansion. Here often a faster convergence of the chiral series is expected. 
Examples of STich parameters are the coefficients Ci of the polynomial C{s,t,u) de- 
fined in Eq. (3.11) after the polynomial ambiguity in the V^{t) is removed by requiring 

^V'(s)=0 (» = 0.1. 2. 3 for/ = 0,2. 

ds" ^ ^ \n = 0,1,2 for /-I. ^ ^ 

These were the quantities used in Ref. [20, 21] to perform the matching of the Roy 
analysis with the SU{2) ChPT constraints. In addition they defined two combina- 
tions of these constants which had very small nonanalytic contributions depending 
on the pion mass when reexpressed with the help of the pion scalar radius. These 
are given by 

Ci = F,'(c2 + 4m^(c3-C4)) , 

F2 

C2 = ^ (-ci + 4mt (C3 - C4)) . (6.9) 
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(a) (b) 

Figure 5: The subthreshold parameters Ci and C2 as a function of the input values L4 
and Lq with the other simultaneously refitted to the K(>4, form factors, (a) Ci calculated 
to order and p^. (b) Idem for C2. Notice that the lowest order value of 1 is outside the 
plotted region. 

They are both equal to one at lowest order. We have shown both Ci and C2 as a 
function of and Lg similarly to the plots shown for and Oq in Fig. 5. Also shown 
are the regions for them obtained in Ref. [21]. 

The plots show again that the preferred values of L4 and Lg are somewhat 
different from zero in order to obtain good agreement. The values of C5 and Cg can 
again be used to check the resonance predictions for two possible combinations of 
the C^. These are in agreement to about a factor of two as expected. The values for 
C3 is in good agreement with the one obtained by Ref. [21] but the contribution 
is of the opposite sign and larger than the contribution. For C4 the agreement is 
marginal, similar to the one for shown in Table 2. 

7. Conclusions 

We have calculated vrvr scattering to NNLO in three flavour ChPT and presented the 
full expressions in App. B. This is the main results of this work. We then presented 
a first numerical analysis and comparison with the low energy constants of order 
as determined from the absolute values of the form factors in Ke4. 

This comparison led to the conclusion that the preferred values of L4 and Lg are 
somewhat different from zero but in size compatible with the error estimates based 
on the arguments from the large number of colours limit done in Ref. [6]. We have 
presented plots of two of the threshold parameters as a function of and Lg with 
the other L[ fitted to the K^a data. The values for the threshold parameters have 
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been presented for fit 10 of Ref. [42] and fits A,B,C of Ref. [32]. The latter are those 
in the region preferred by the various scalar form factor constraints. The convergence 
of the three flavour ChPT series for tttt scattering seems reasonable and is similar to 
the one for the two flavour case after the masses and decay constants have been put 
to their physical values as was done here. 
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A. Integrals 

The one loop integrals we use are well known. The expressions for them can be 
found in [33] and [36]. The two loop integrals of sunset type were derived in [33] 
using the results of [51] for the subtraction constants. They are denoted by in 
the expressions. 

The two loop integrals of vertex type are evaluated using the methods of [52]. 
Their precise definition and method of evaluation can be found in [36]. They are 
denoted by V in the expressions. 

B. Analytical results 

We now write the amplitude in the form (3.9) where we expand all functions in the 
ChPT expansion which we label by a superscript (n) denoting the order in the 
expansion. The FORM files with this output will be made available in [53] . 
The lowest order result is [1] 

V'^'\s)^0, C{s,t,uf^^-^{s-ml) . (B.l) 

The order p'^ result agrees with the one shown in [10] and with the one in [9] up to 
a few misprints there. The result is 

327rF,^ (s) = +B{ml, ml, s) (-2 ml s + 1/2 + 2 s^) 

+B{ml,ml,s) {3/8s^)+B{ml,ml,s) {l/6ml) , 

327rF^V'^^\s) ^ +B(ml,ml,s) (-2/9 + 1/18 s) 

+B{ml,ml,s) (-l/9m^ + l/36s) , 

32ttF^V^^''\s) = +B{ml,ml,s) (-2 s + 1/2 + 2 s^) 

+S(m^, ml, s) (3/8 s^) + B{ml, mj, s) (l/6 m^^) , (B.2) 



16 



and for the polynomial part 

C{s, t, u)^""^ = (2/3 ml m\ - 1/3 ml s + 2/3 m\ - 1/2 m\ s - 1/8 + 1/24 S^) 

+A{ml) (-4/3 ml + s)+ A{m]^) (-2/3 ml + 1/2 s) 
-32 ml L\s - 16 L^s + 16 + 8 L^s 

+32 + 16 ml - 32 m^ LI - 16 m^ L; 

+32 + 16 + 8 L[s^ + 2 L^s^ + 2 L^^^ + 4 L^s^ (B.3) 

Here we used the notation 5 — t — u to have a shorter expression. 

The order expressions are significantly longer. We split them in several parts 

F^C^''\s,t,u) = Ccis,t,u) + Cj^{s,t,u) + Ceis,t,u) , 

32nF^ V'^'\s) = Vlis) + V^{s) + V^{s) . (B.4) 

The C[ only contribute to the polynomial part. 

Cc{s,t,u) = +s5^ {+6C; + 2Cl) 

+s ( - 64 ml m\ Q + 64 ml m\ - 256 ml m\ Cl^ + 32 ml m\ 
-96 ml CI - 192 ml CI + 96 ml CI + 32 ml CI - 32 ml CI - 32 ml CI 
-64 ml Clj - 32 Cg" + 32 C[o + 32 m^ C{^ - 64 C[^ 
+32 C[4 + 48 + 64 + 32 C{.j - 32 C2'5 

+32 ml Cl^ + 64 ml C^, - 64 m^^ C^, - 32 Cg^o) 

( + 48 C[ + 96 ml C\ - 48 CI - 40 Q + 8 
+8 + 16 ml CI, + 8 Cg" - 4 ml C[o - 4 C[i + 8 C\^ 

-8 - 8 ml C^, - 8 ml C^, + 20 CI,, 

+ 16 Q - 8 m\ C[i + 48 m\ ^[3) 
+5^ (-8 C[ - 16 C^' + 10 CI + 14 CI) 

+8^[-SmlCl + Aml C{, + 4 ml CI, - 8 C{, - 8 ^[3 

+8 ml CI, + 8 ml C^^ - 4 + 8 m^, C[i - 16 m^, ^[3) 

+64 ml ml CI - 64 ml m\ C{, + 192 ml m\ C{^ - 64 C{, 

+64 C2^o + 384 ml m\ C^, + 64 ml m\ CI, + 64 ml CI 

+128 ml CI - 64 C^ + 32 + 32 ml CI + 64 ml Clj + 32 m^ CI 

-32 C[o - 32 ml C{, + 32 ml CI, - 32 ml Cl^ - 64 C{^ - 96 

-128 ml C[6 - 64 ml C^^ + 96 C[^ + 160 C2^o + 192 ml C'„ 

-64 ml CI, - 128 C^^g + 64 ml C^, + 96 ml C^, + 160 ml CI, . (B.5) 

Ci^{s,t,u) = +— ^ ( + 256/27 m^m^L^s + 1048/81 m^m^L^s- 8 m^m^L^s 



17 



-92/9 L\s^ + 28/9 LIS"^ - 386/27 mlLls^ -2/9 ml 1^6^ 
-320/81 ml Lls^ + 4/3 L16^ - 8/3 + 8/9 ml L\5^ -2ml ^5^' 
+2/3 ml L16^ + 128/27 ml m\ - 736/81 m^ m\ LI + 32/3 m^ m% 
+176/3 ml Lis + 808/27m^ L^s + 2008/81 ml L^s - 160/9 ml L^s 
- 16/3 m^L^s- 448/9 m^L^ - 224/27 m^L^ - 1760/81 m^L^ 
+256/9 ml LI + 32/3 ml LI - 8/3 m\ L^s^ - 8/3 m| L^s^ + 2/3 L^f^^ 
+1/9 L^s 5^ + 1/3 L[s^ + 5/18 L^s 5^ + 13/6 L^s^ + 1/9 L^s^) 

+A(ml) ( - 16/3 ml L\s - 16/9 ml L^s - 40/27 ml L^s + 32/3 ml L\ 
+16/9 LI + 64/27 - 16/3 ml LI - 8/9 ml L; + 32/3 ml 

+ 16/3 m^L^) 

+A(m2 ) ( - 32/3 ml L\s + 184/3 ml L^s - 16/3 ml Lis + 64 L^s 
+32 + 64 ml L\ - 32 ml L^ + 32 ml L^ - 464/3 ml LI 

-232/3 ml LI + 160 m,l L^ + 80 L^ - 12 L[s^ + 4 L^5^ 
-50/3 L^s^ + 38/3 L'^S^ - 6 L^s^ + 2 L^5^) 

+A(m^) ( - 32 L^s + 16/3 L^s - 52/3 ml L^s + 16 ml L^s 

+16 Lis + 64 + 80/3 ml L^ - 64 ml L^ - 104/3 ml LI 

+64 + 32 ml L; - 4/3 L^s' + 4 L^^' + 5/3 L^s' + 1/3 L^S^) 

+128 L^L^s - 512 ml m\ L^L^s + 256 ml m\ {L^f s 

-256 ml ml L^L^s - 256 ml m\ L\Ll + 2048 ml m^ L^L^ 

+512 ml ml L^L; - 512 ml m\ {Llf + 512 ml m^ LIL^ 

-1024 ml ml L^L; - 2048 ml ml {Llf + 192 ml L^U^s 

-256 L^L^s - 256 ml L^L^s + 128 ml {L^f s - 128 ml L^L^s 

-128 LlLls + 64 (L^)' s - 384 L^L^: + 1024 ml L\Ll 

+768 ml LIL; - 256 {Llf + 768 L'^L;: + 512 ml LlLl 

-128 ml {Llf- 1536 L^L^ - 1024 ml {L^f - 512 {L^f . (B.6) 



C6(s,i,ii) = +-^ A(m2) ( + l/36m^s - 1/9 m^) 



^(m^) ( + 8/9 - 89/36 s + 29/18 - 2/3 s 



167r2 
1 

^16^ 
+5/24 s'- 1/72 5') 
1 



A(m|^) + 4/9 ml m^ - 37/36 mls + 3/4 - 1/3 m^ s 



167r2 

+1/245^+1/72 (^2) 
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+ (— ^) (-83/216 mj^pi^ s — 1375/432 m|- s — 2/9 m^^ 

+13/32 TT^ - 19/288 pi^ S'^ + 1751 /432 ml - 289/432 5^ 
+103/648 tt^ + 167/432 - 337/216 tt^ s 

-10933/864 s + 349/324 tt^ + 6013/864 ml + 1/8 tt^ 
-1/144 TT^ 5^ + 13/12 - 7/72 m\ 5^ + 1/6 s 

-1/96 s 5^ - 11/96 s^) 

+A{mlf ( + 1/12 ml m'^ s - 2/9 m^ m'^) 

+A{ml) A(m^) ( - 8/9 ml + 2/3 s) 

+A{ml)^ ( + 39/16 m;' - 19/48 m;' 5' + 11/2 ml - 26/3 s) 
+]4(m^)2 ( - 29/12 ml m'^ s - 2/9 ml + 5/4 m^ m^^ + 3/4 m^ 
-l/2AmK^5'^ + l/Qs) 

+H^\ml, ml, ml, ) ( + 5/2 m^s- 5/2 m^) 
+H^'(ml, m%, m\, ml) ( - 15/8 s + 15/8 m^) 
+H^'iml, ml,ml ml) ( + 1/6 s - 1/6 m^) 

, ttFi f 2 2 2 2\ ( , o 2 2 o 4 2\ 

+// [mj^,m^,mj^,mj + Sm^mj^ s - Sm^mj^j 
+H^'{ml, ml, m\, ml) ( - 5/2 s + 5/2 m^) 
+H^'{mlj,mj(,m%,ml) ( + 3/2m^m|-s — 3/2m^m|^ + 13/8 m^s 
-13/8 m^) 

+H['{ml, m\, m^, m^) ( + 3 s - 3 m^) 

+H[\ml, m\, m%, ) ( + 2 s - 2 m^) 

+H[\ml, m\, ml,ml) (^-2mls + 2 m^) 

+H^^{ml, ml,ml,ml)( + 9mls-9 ml) 

+H^;{ml, m\, ml, ml) ( - 9/8 s + 9/8 ml) 

+iff/(m^, ml, m\, ml) ( + 9 s - 9 ml) 

+i/f/(m;, ml, m\, ml) ( + 27/8 ml s - 27/8 ml) 

+H'"{ml, ml, ml, ml) [ + Amis- 17/6 - + 1/3 S^) 

+H^{ml, m\, m\, ml) ( - 133/36 s + 10/3 m^ - 5/16 + 7/16 5^) 

+i/^(m^,m2,mj,m2) ( + l/18m^) 

+H^{ml, ml, ml, m^) ( - 1 /12 ml m\ + 53/18 s - 4/3 m^ 
+13/72 m|s- 1/4 s^- 1/4 5^) 
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+H^{ml,ml,ml,ml) (^ + l/12mlm% - 13/72 m^s) 
+H^{m\, m\, ml, ml) ( - 7/18 m^ + 3/16 - 3/16 5^) 
+H^{ml, mj,, ml, ml) ( - 11/8 ml s + 119/36 m^ - 45/64 + 27/64 5^) 
+H[{ml, ml, m\, ml) ( + 7/4 s - 3/2 m^ - 3/16 + 1/16 5^) 
+i/f (mj, m^, m\, ml) ( + 11/8 s - 15/4 + 45/64 - 27/64 
+i/fi(m2 , ml, ml, m^) ( + 5 s - - 3/2 + 1/2 5^) 
+H[i{ml, ml, ml, ml) ( - 13/8 mls + 1/8 m^ + 15/16 - 5/16 (5^) 
+//2^i(m^, ml, ml, ml) ( + 253/24 mis- 45/4 + 3/8 - 1/8 5^) 
+H^, {ml, ml, ml, ) ( - 37/24 ml s + 17/4 - 3/2 + 1/2 5^) 
+H^^{ml, ml, ml, ml) ( + 3/2 ml s - 9/2 rr4 + 45/32 - 15/32 5^) 
+H[^{ml,ml,ml,ml) (^ + 21/8mls - 3/8ml-9/8s^ + 3/85^) . (B.7) 

Vl°(s) = +B{ml,ml,s) ( - 728/3 m^L^s^ - 376/3 m^L^s^ - 364/3 m^L^s^ 

+64 ml Ll-s^ + 32 ml Lls^ + 1024/3 m^ L[s + 608/3 m^ L^s + 512/3 m^ L^s 
-224 L^s - 112 L^s + 320 L^s + 160 L^s - 352/3 ml L[ 
-224/3 ml LI - 176/3 + 96 + 48 ml LI - 160 ml L^ 

-80 LI + 176/3 L^s^ + 112/3 U^s^ + 88/3 L^s^) 

+B(ml, ml, s) (^ + 96 ml ml L\s + 32 ml ml L2S + 32 ml ml L^s 
-96 ml ml L^s - 24 ml ml Lis + 96 ml ml L^s + 48 ml ml L^s 
-48 ml Lls^ -Sml L^s'^ - 14 ml L^s"^ + 24 ml L^s"^ + 6 L^s^ 
-48 L^s^ - 8 L^s^ - 14 L^s^ + 24 ml L^s'^ + 24 L\s^ 
+8 L'^s^ + 8 L^s^) 

+B{ml, ml, s) ( - 64/3 ml ml L\s - 32/9 ml ml L^s - 128/27 ml ml Lis 

+32/3 ml ml L^s + 8ml L[s^ + 8/3 ml L^s^ + 20/9 ml L^s^ 

+128/3 ml ml L[ + 128/9 m^ ml L^ + 320/27 m^ ml LI - 128/3 m^ ml LI 

-64/9 ml LI + 128/3 ml L^ - 128/3 ml - 32/3 L^s 

-16/9 L^s - 64/27 ml L^s + 16/3 ml Lis - 32/3 L^ - 32/9 

-80/27 ml LI + 32/3 ml LI + 16/9 - 32/3 L^: + 128/3 ml 

+16 m^L^). (B.8) 

V^{s) = +B{ml, ml, s) ( + 85/24 ml s" - 25/3 s + 10/3 ml - 5/12 s^) 
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+B{ml, ml,sf (^-3mls^ + 3/2mls-l/Aml + 2. 
+B{ml, ml, s) B{ml, m|, s) ( - 3/8 ml + 3/4 s^) 
+S(m^, ml, s) B{ml, mj, s) ( + 1/3 s - 1/6 m^) 
+S(m^, m^, s) ( - 20 s ]4(m^) - 10 s A(m^) + 8 ]4(m^) 

+4 ml A{ml) + 8s^ A{ml) + 4 A(m^) 



+S(m^, m^, s) ^ ( - 2 s + 1/2 ml + 1/2 m\ - 1/8 s^) 

+5(m^,m^,s)'( + 9/32s=^) 

+B{m\, m\, s) B{ml, mj, s) ( 
-1/3 ,s + 3/8 ml s^) 

+S(m^, m^, s) ( - 3/4 ml ]4(mJ) s - 9/4 s ]4(m2 ) - 3/2 s A{mj^) 
+3/2 Z(m^) + 3/4 A{mj,)) 

+B{ml, ml, s) ^ ( + 2/9 mlm^s- 1/24 ml - 8/9 m^ ml + 1/9 m^ s 
+2/9 ml) 

+B(ml, ml, s)' ( + 4/27 - 7/108 ml) 
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+B{ml, ml, s) ( + 2/3 m^ A{ml) - 2/3 m^ A{m'^] 

+B\ml,ml,s) — ^ ( - lOmlm^s + 59/4 s'^ + 4mt 

-335/9 s + 136/9 + 4 - 3/2 s^) 

+B'^{m\,m\,s) — ^ ( - Qmlm\s + 15/8 m^s^ - 2m^s + 9/4 
167r ^ 

-9/16 s^) 

+F(mJ, mj, s) ^ ( + 1/3 mlm\s- 1/4 ml - 4/3 m\ + 2/3 s 
-2/9 m^). (B.9) 



V^{s) = +Vo(m^,m^,m^,m^,m^,s,m^) ( + 4m^ - 4m^ s + m^) 

+Vb(m^, m^, m^, m^, ml, s, ml) ( + 1/2 ml - 9/4 s + m^) 
+yo(m^, m^, mj, ml, ml, s, ml) ( + 2/9 s - 1/9 m^) 
+\/o("^i, "^i, ml, s, ml) ( + 9/2 s) 

+Vo{m\, m\, m\, ml, ml, s, ml) ( - 3/8 + 2 m^ s) 
+Vb(mJ, ml, ml, ml, ml, s, ml) ( + 2/9 m^) 
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+Vo(m^^, m^^, m^, m^, m^, s, ml) + 1/4 rr4 s + 1/3 m^^ 

iTA/'2 2 2 2 2 2\/ 22iri/o4 o6\ 

+Vu (m^, m^,mK, tuk, m^, s, mjl^-m^s + 9/2 m^) 

iTA/2 2 2 2 2 2\/'io/o22 c4N 

+\/ii(mJ,mJ,m^,m^,m^,s,m^) (+ l/2m;i6- - 2m^) 
+Vis{mj^,m'j^,ml,m^i^,ml, s,ml) + 21/4 - 9m^ s - 3/4 s^) 
+Vi3{mj^,mj(,mj(,m^,rnl,s,rnl) 15/4 m^s^ - 6m^s - 9/16 s^) 
+yi4(m|-,m|^,m^,m^,m^,s,m^) ( + 6m^s^ - 3/2 s^) 

+Vu{ml,m^j^,m^j^,ml,ml,s,ml) ( + 9/2 - 9/8 s^) 

iT//2 2 2 2 2 2\/ /|22iio4 o6\ 

+ V22{'m^,m^,m^,m^,m^,s,mJ y - Am^s + 18 m^s - 8m^j 
+V'22(m^,m^,m^,m^,m^,s,m^) (+ l/2m^s^ - 9/4 m^s + m^) 
+^22(m^, mj^, ml, m^, s, m^) ( - 3 ml + 12 s) 
+V22{m\, m^, m|^, m^, ml, s, ml) ( - 9/8 ml + 9/2 s) 
+V22{ml, m\, m\, m\, ml, s, ml) ( - 3/4 s + 3 m^) 
+V23{ml,ml,ml,ml,ml,s,ml) (- 18 + 8 s + 4 s^) 
+V'23("^^, ml;, ml;, ml, s, ml) ( + 9/4 - s - 1/2 s^) 
+V23{ml,mK,ml,ml,ml,s,ml) ( - 12 + 3 s^) 
+V^23(rni, "^K, "^K, 5, "^^) ( - 9/2 + g/g 

iTA/2 2 2 2 2 2\/|0//f22 o4\ 

+V23{m^,m^,mK,mj^,m^,s,mJ + 3/4 m^s -3m^sj 
+V'26("^^,^«|-,?«^"^x,"^^s,"^^) ( - 6m^s^ + 12m^s + 3/4 s^) 
+V2e{m1,m1,m1,ml^,ml, s,ml) — 9/2 ml s"^ + 9m^s + 9/16 s^) 
+V27(ml, ml, ml, ml, ml, s, ml) ( - 12 + 3 s^) 

+V27{ml, ml, ml, ml,ml,s, ml) ( - 9m^s^ + 9/4 s^) 

IT//2 2 2 2 2 2 \ f 22,0 /0 3\ 

+V29[mj^,mj^,m^,mj^,m^,s,m^)y-Qm^s +3/2s j 

+1^29 ("^/sr> "^is:> -s, "^^) ( - 9/2 ml + 9/8 s^) 

+V2io(ml,ml,ml,ml,ml,s,ml) ( + 40m^ s - 16m^ - 16 s^) 
+V'2io(?7i^,?7i^,?7i/s:,"^^,?7i^,s, "^^) ( + 20m^s - 8m^ - 8s^) 

+V2io{ml, ml, ml, ml, ml, s, ml) ( + 9/2 s - 9/4 s^) 
+V2io{ml, ml, ml, mj, ml, s, ml) ( + 9/2 s - 9/4 s^) 
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+V2ni'ml,ml,ml,ml,ml, s,ml) — 16 m^s^ + 40m^s — 16 m^) 

iTA/S 2 2 2 2 2\/ o 22, on 4 o6\ 

+^211 (m^, m\,m\,m\,m\, s, ml) ( - 9/4 + 9/2 s) 
^A^2ii{m\,m]^,m\,m'l^,ml,s,ml) ( - 9/4 + 9/2 s) 

IT//2 2 2 2 2 2\/|^22 o3\ 

+V2i2{m^,m^,m^,m^,m^,s,mJ {^ + Am^s -8s j 
+V2i2{ml, ml, m\, m\, ml, s, ml) ( + 2 - 4 s^) 
+1412 (m^, m\,ml,m\,ml, s, ml) ( - 9/8 s^) 
+^2i2(m^, m%, m\, mj, m^, s, ml) ( - 9/8 s^) 
+V2i3{ml,ml,ml,ml,ml,s,ml) ( + 4m^s^ -8s^) 

IT//2 2 2 2 2 2\/io22 ,i3\ 

+V2i3{m^,m^,mj^,mj^,m^,s,mJ {^ + 2m^s -4s J 
+V^2i3(m^, m\,ml,m\,ml, s, ml) ( - 9/8 s^) 

+14i3(m^, ml, m\, m% ml, s, ml) ( - 9/8 s") . (B.IO) 

Vl{s) = +:B(m2 , ml, s) ( + 32/9 ml L\s - 16/9 L^s + 16/9 L^s + 16/9 ml L^s 

+8/9 ml Lis- 64/9 - 32/9 m^ - 8/9 L^s^ + 4/9 U^s^ - 4/9 L^s^) 

+B{ml, ml, s) ( - 16/9 m^ L; + 4/9 rr^ L^s + 8/9 L^s - 2/9 L^s^) . 

(B.ll) 

Vi{s) = +S(m^, m^, ( + 4/27 - 11/27 s + 28/27 - 1/27 s 

+1/27 s^) 

+B{ml,ml,s)^ [-2/27 mis + A/27 ml + 1/108 s^) 

+B{ml,ml, s)B{ml,ml, s) ( + 4/27 m^m^ - 1/27 m^s - 1/27 m^s 

+1/108 s^) 

+B{ml,ml, s) [ - 8/27 ml A{ml) - A/27 ml A{ml) + 2/27 sA{ml) 
+ l/27sA{ml)) 

+B{ml, ml, s) ( + 10/27 ml ml - 5/54 s - 5/54 mls + 2/27 ml 
+1/54 s'') 

+B{ml, ml, s)^ [ - l/5Aml s + 1/27 ml + 1 /432 s^) 

+ B {ml, ml, s)[- 4/27 ml A{ml)- 2/27 ml A{ml) + 1/27 sA{ml) 
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+ 1/54:8 A{m%)) 

+B'{ml, ml,s)^{- 7/18 ml s + 10/9 + 1/36 s') 

+B\m\, ml, s) ( + 1/3 ml ml - 1/12 ml s - 1/18 ml s + 1/72 s^) . 

(B.12) 



V^{s) = +Vo{ml,ml,ml,ml,ml,s,ml) (- l/3m^s) 

, , ml, ml, ml, s, ml) ( - 2/3 s + 1/6 s^) 
+Vo{ml,ml,ml,m'^,ml, s,ml) (- l/3m^s + l/8s^) 

+1^11(772^, m^,m^,m^,m^,s,m^) ( + 1/3 m^s + 1/3 m^) 
+Vii(m^,m^,mJ,mJ,m^,s,rr4) (+ 1/27 m^) 
+^ii(m|^,m^,m^,m^,m^,s,m^) ( + 2/3 m^s + m^ - 1/4 s^) 
+Vn(ml,ml,ml,m'^,ml, s,ml) ( + l/3m^ s + 4/9m^ - 3/16 s^) 
+Vi3{ml,ml,ml,ml,ml,s,ml) ( + 2/3 m^s- 1/6 s^) 
+Vi3{ml,ml,ml,m'^i,ml,s,ml) ( + 1/2 m^s - 1/8 s^) 
+Vu{ml, ml, ml, ml, ml, s, ml) ( - s + 1/3 s^) 
+Vi4{ml,ml,ml,m'^,ml, s,ml) ( - 2/3 s + 1/4 s^) 
+V22(ml, ml,ml,ml,ml,s,ml) - A/3 m^) 
+V22{ml, ml, ml, ml, ml, s, ml) ( + 2/3 s - 8/3 m^) 
+1^22 {ml, ml , m,l , mf^ , m,l ■s,ml) + 1/2 s — 4/3 m^) 
+V23{ml, ml, ml, ml, ml, s, ml) + A/2> ml s) 
+V23(ml,ml,ml,ml,ml,s,ml) ( + 8/3m^s - 2/3 s^) 
+V23{ml,ml,ml,m'^^,ml, s,ml) ( + 4/3m^s - 1/2 s^) 
+V2Q{ml,ml,ml,ml,ml, s,ml) ( - 1/6 s - 2 + 1/6 s^) 
+V2&{ml,ml,ml,m1^,ml,s,ml) ( - 1/6 m^s -4/3m^ + 1/8 s^) 
+V27(ml, ml, ml, ml, ml, s, ml) ( + 2 s - 2/3 s^) 
+V27(ml,ml,ml,m'^,ml,s,ml) (^ + 4/3m^s - 1/2 s^) 
+V2g{ml,ml,ml,ml,ml,s,ml) (+ 11/3 m^s - s^) 
+V2g{ml,ml,ml,ml,ml,s,ml) ( + 8/3 m^s - 3/4 s^) 

, -tr / 2 2 2 2 2 2\/|/i2 \ 

+V3i(m^, m;^,m^,m;^,m^,s,m^) (^ + 4m^ - sj 
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+V3i{mj^,mj^,m]^,ml,ml,s,ml) ( + 3m^ - 3/4s) 
+V33{ml, ml, ml, ml, s, ml) ( - 2/3 s + 4/3 m^ + 1/12 s^) 
+V33(ml, ml, ml, mj, ml, s, ml) ( - 1/2 s + + 1/16 s^) 
+V2i3iml,ml,ml,ml,ml,s,ml) ( + 8/3m^s - 2/3 s^) 

+V2n{ml,ml,ml,ml,ml,s,ml) ( + 4/3m^s - 1/3 s^) 
+V2i3{ml,ml,ml,ml,ml,s,ml) ( + s - 1/4 s^) 

I / 2 2 2 2 2 2\ ( , 2 i //( 2\ 

+ V/2i3("^^, "i;^, m^, rn^, s, mj + s - 1/4 s j 
+y3ii(m|-, , m^, m|^, ml, s, ml) ( - 2/3 s + 8/3 m^) 

+V3ii{ml,ml,ml,m'^^,ml,s,ml) (- l/2m^s + 2m^) 
+V3i3(m^,m^,m^,m^,m^, s,m^) ( - 8/3ml s + 2/3 s^) 
+V3i3{ml,ml,ml,m'^,ml, s,ml) ( - 2 s + 1/2 s^) 
+y3i4(m|^, m|^, ml, ml, ml, s, ml) ( - 8/3 s + 2/3 s^) 
+y3i4(m|^, m|-, ml, m^, ml, s, ml) ( - 2 s + 1/2 s^) 
+^3i6("^i,"^i,"^^"^x,"^^,s,m^) (- 16/3 m^s + 4/3 s^) 

iT^/2 2 2 2 2 2\ ( A 2 i2\ 

+V3ie{mj^,mj^,mj^,m^,m^,s,mJ \^-Am^s + s j 
+V3i7{ml,ml,ml,ml,ml,s,ml) ( + 8/3 - 4/3 s) 
+V3i7{ml, ml, ml, ml, ml, s, ml) ( + 4/3 ml - 2/3 s) 
+V3i7{ml, ml, ml, ml, ml, s, ml) ( + - 1/2 s) 
+V3i7{ml, ml, ml, m\, ml, s, ml) ( + - 1/2 s) 
+1^321 {ml, ml, ml, ml, ml, s, ml) ( - 4/3 s + 8/3 m^) 
+V32i(ml,ml,ml,ml,ml, s,ml) ( - 2/3m^ s + 4/3m^) 

, 17- / 2 2 2 2 2 2\ f 1 /o 2 I A\ 

+V321 [nriK, m^, m^,mK,m^, s, mj - 1/2 s + m^j 

+V32i{ml,ml,ml,ml,ml,s,ml) (- 1/2 m^s + m^) 
+V'323("^^, ml, ml, ml, ml, s, ml) ( - 2/3 s^) 

, -tr / 2 2 2 2 2 2\ f i/o2\ 

+V323{m^:m^,mK:mK,m^,s,mJ l/3s j 
+V323(ml, ml, ml, ml, ml, s, ml) ( - 1/4 s^) 
+V323iml, ml, ml, mj, ml, s, ml) ( - 1/4 s^) 

+V325{ml,ml,ml,ml,ml,s,ml) ( - 8/3 s) 
+'^325("^^, m^, m^, m^, m^, s, m^) ( - 4/3 s) 
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2 2 2 2 2\ ( 2 

+ V325("^A-, ^^^X, ^^i^, ^^K^ "^TT^ "^Tt) ~ "^TT 

iTA/'2 2 2 2 2 2\ ( 2 \ 

+Vz2b[mK, m^, rriK, m^, m^, s, m^) - sj 
+V326(ml,ml,ml,ml,ml,s,ml) (- 16/3 s + 4/3 s^) 

+V^j26(m^, m^, m^, m%, ml, s, ml) ( - 8/3 s + 2/3 s^) 

+V'326("^i, "^i, ml, s, ml) [-2ml + 1/2 s^) 

+'^326("^^> ^% s,ml) [-2mls + 1/2 s^) . 



(B.13) 



(s) = +B{ml, ml, s) ( + 64/3 ml L\s^ + 176/3 ml L^s^ + 32/3 ml L^s^ 

+ 16 L:js2 + 8ml L^s^ - 128/3 L[s - 352/3 L^s - 64/3 m^ L^s 
-32 ml Lis - 16 Lis - 64 Lgs - 32 ml L^s + 128/3 ml L\ 
+256/3 + 64/3m^ L^ + 128m^ Ll + UmlLl - 16/3 L^s^ 
-32/3 L^s^ - 8/3 Lls^) . (B.14) 



Vi (s) = +B(m2 , , s) ^ ( - 5/6 + 8/3 m^ s - 8/3 m^ + 1/12 s^) 
+S(m^, ml, sf ( + 3/2 - 3 s + 2 - 1/4 s^) 
+B{ml, ml,s)(-8mls A{ml) -Amis A{ml) + 8 m^^ A{ml) 

+4 ml A{m\) + 2 A{ml) + Z(m^)) 

+S'(m^,m^,s) [ - Amlm\s - Ami s"^ + Am^mj^ + 73/9 m^ s 
-56/9 m^ + m^s^ + 3/4 s^) . (B.15) 



V^(s) = +Vo(ml,ml,ml,ml,ml,s,ml) ( + m^s^ - 7m^s + lOm^) 

+yo(rn^,"^^,"^i,?74,?7i^,s,"^^) ( + 1/2 m^s^ - 3m^s + 4m^) 
+Vo{ml,ml,ml,ml,ml,s,ml) (- 1/9 s + 2/9 m^) 
+l^ii(m^,m^,m^,m^,m^, s,m^) ( - 3 + 18 s - 24 m^) 
+Vii{ml, ml, m\, m]^, ml, s, ml) ( - + 6 s - 8 m^) 

^-V22{ml,ml,ml,ml,ml,s,ml) ( + 2m^s^ - 12m^s + 16m^) 
+V22(m^, m^,m^,m^,m^,s,rr4) ( + 1/2 m^s^ - 3m^s + 4m^) 
+V2'i{ml,ml,ml,ml,ml,s,ml) (+ 12m^s^ - 16m^s - 2s^) 
^-V2z{ml,ml,m\,m\,ml,s,ml) ( + 3m^ - 4m^ s - 1/2 s^) 
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+V2io{ml, ml, ml, ml, ml, s, ml) ( + 16 s - 16 - 4 s^) 

, ir / 2 2 2 2 2 2\/io2 o4 o2\ 

+V2ii(ml, ml, ml, ml, ml, s, m^) ( - 4 + 16 s - 16 mfj 

, -fr / 2 2 2 2 2 2\ f o22io4 o6^ 

+V2ii(m^, m^, mj^, m^^, m^, s, mj y-2m^s +8m^s-8 m^j 

IT//2 2 2 2 2 2\/|/|22 o3\ 

+ V2i2(m^, m^, s, mj i^ + Am^s -2s j 

, ir / 2 2 2 2 2 2\/io22 3\ 

+ V212 ("^^, A- >"^/^'"^7r)S,m^) (^ + 2 m^s -s j 

iT^/'2 2 2 2 2 2\/ij22 o3\ 

+1^2i3(m^,m^,m^,m^,m^,s,m^) (^ + 4m^s -2s j 

+V^213("^^, ml, m\, m\, ml, s, ml) ( + 2 - s^) . (B.16) 
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